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Deformation  Criteria  for  the  Direct 
Manipulation  of  Free  Form  Surfaces 


Stephane  Guillet  and  Jean-Claude  Leon 


Abstract.  The  approach  proposed  here  is  based  on  coupling  a mechan- 
ical model  to  the  input  surface  geometry  provided  by  the  designer.  The 
mechanical  model  is  based  on  a set  of  bar  networks  subjected  to  tension 
forces  only.  Constraints  are  specified  by  the  designer  to  express  the  defor- 
mation behaviour  of  the  surface  in  the  area  of  interest  and  define  functional 
dimensioning  objectives  to  perform  a direct  manipulation  of  the  surface. 
Generally,  this  process  leads  to  a globally  underdetermined  system  of  equa- 
tions, i.e.  the  number  of  unknowns  (external  forces)  is  significantly  greater 
than  the  number  of  equations  generated  by  the  designer’s  constraints.  To 
this  end,  a minimization  problem  is  formulated  which  expresses  various  de- 
formation behaviours.  In  contrast  to  difference  of  deformation  approaches 
based  on  mechanical  models  like  membrane  models,  finite  element  models, 
which  solely  rely  on  strain  energy  minimization  criteria,  the  approach  pro- 
posed here  provides  the  designer  various  criteria  to  help  him/her  create 
different  deformation  behaviours  like  an  area  minimizing  criterion,  ex- 
pressing a minimum  change  of  the  shape  in  the  deformed  area,  expressing 
a deformation  with  slowly  varying  curvature  in  the  deformed  area,  pro- 
viding a deformation  behaviour  which  allows  to  approximately  preserve 
the  section  of  pipe-like  surfaces  subjected  to  bending  deformations.  As 
depicted,  multiple  criteria  help  the  designer  express  various  deformation 
behaviours  which  are  required  during  a design  process. 


§1.  Context  of  Surface  Deformations 

The  shape  modification  of  an  object  during  a design  process  depends  on  the 
context  of  this  process.  In  the  field  of  mechanical  engineering  design,  con- 
straint requirements  can  be  either  aesthetic  [3,7]  or  functional  [1], 

Without  adequate  3D  modification  tools,  the  surface  deformation  leads 
the  designer  to  tedious  manipulations,  i.e.,  displacements  of  numerous  control 
polyhedron  vertices,  chain  modifications  of  patches  or  surfaces,  etc. 

The  basic  aim  of  these  deformation  tools  would  be  to  provide  the  user  an 
easy  and  intuitive  control  of  the  surface  shape.  Their  parameters  should  be 
automatically  related  to  the  parameters  governing  the  deformation  process. 
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Such  approaches  fit  into  a class  [3,5,9]  that  helps  a designer  shape  the  overall 
object,  but  they  are  not  suitable  for  generating  free- form  surfaces  which  accu- 
rately match  geometric  constraints  involving  functional  parameters.  A second 
class  of  3D  modification  tools  covers  the  approach  of  Celniker  [2],  Kondo  [6], 
Light  [8]  and  Welch  [10].  These  types  of  tools  fall  into  the  domain  of  paramet- 
ric or  variational  design  tools.  Among  these  approaches,  some  [6,8]  focus  on 
parametric  or  variational  models  applied  to  2D  models.  Others  [2,10]  perform 
a surface  deformation  subjected  to  constraints  such  as  prescribed  curvature 
or  surface  rectitude  using  a membrane  model  which  cannot  provide  some  de- 
formation modes  like  bending. 

The  approach  presented  here  fits  into  this  last  category.  Similar  to  the 
approach  of  Celniker  and  Welch,  the  current  one  also  uses  a mechanical  model. 
However,  its  formulation  is  simpler  than  Celniker  and  Welch’s,  and  thus  it  is 
easier  to  manipulate  and  it  allows  generation  of  isotropic  and  anisotropic  de- 
formations. The  approach  introduced  here  is  a new  development  around  a 
free- form  deformation  method  [4],  The  scope  of  the  present  work  focuses  on 
the  introduction  of  a set  of  deformation  criteria  which  cannot  be  provided  by 
membrane  models  or  other  mechanical  models  subjected  to  small  displace- 
ments and  linear  behaviour  material  law  hypotheses. 

§2.  Principle  of  the  Parametric  Deformation 

Before  studying  the  deformation  criteria  presented  here,  it  is  suitable  to  sum- 
marize the  objectives  and  the  constraints  related  to  the  parametric  deforma- 
tion process.  The  features  of  the  current  work  are  the  following: 

• the  treatment  of  configurations  involving  multiple  trimmed  free-form  sur- 
faces based  on  a B-Spline  model, 

• the  direct  manipulation  of  the  geometry  through  a small  number  of  pa- 
rameters to  allow  an  easy  and  intuitive  control  of  the  surface  shape, 

• the  possibility  for  the  user  to  create  local  or  global  deformations  of  the 
geometry  and  to  obtain  different  solutions  with  one  set  of  geometric  con- 
straints, 

• the  fast  computation  which  allows  an  easy  integration  of  the  parametric 
deformation  tool  into  an  integrated  design  process. 

In  the  context  of  the  approach  introduced  here,  different  constraints  reduce 
the  complexity  of  the  problem: 

• the  surface  patch  decomposition  is  preserved,  i.e.,  degrees,  nodal  se- 
quences and  topology  are  kept  constant, 

• C°  continuity  between  patches  is  maintained.  G°  and  G 1 continuities 
are  approximated  along  the  trimming  lines.  To  this  end,  a discretization 
process  is  applied  to  these  lines  without  modification  of  their  degree, 

• trimming  lines  on  the  surface  are  kept  unchanged  into  their  parametric 
space.  A trimming  line  is  defined  as  a set  of  connected  trimming  curves. 

The  aim  of  the  parametric  approach  is  to  deform  a set  of  trimmed  free-form 
surfaces  subject  to  geometric  constraints.  The  geometric  constraints  are  cur- 
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New  external  forces  Mechanical  features 


Fig.  1.  Main  steps  of  the  parametric  modelling  process. 

rently  restricted  to  the  control  of  the  dimensions  of  an  object  though  other 
categories  of  constraints  can  be  set  up  to  achieve  other  functions. 

Prior  to  the  description  of  the  minimization  criteria  used  here,  it  is  suit- 
able to  describe  how  geometric  and  mechanical  features  fit  together.  The 
resulting  surface  geometry  is  obtained  from  the  initial  one  through  five  steps 
(Figure  1): 

• step  one  is  devoted  to  the  creation  of  the  geometric  features.  These 
features  help  specify  some  target  parameters  of  the  surface  shape,  the 
deformed  area  and  the  continuity  conditions  between  trimmed  patches. 
In  the  case  of  Figure  1,  one  feature  is  generated  by  the  user.  Currently, 
the  design  constraints  can  be  displacements  of  points,  lines;  tangency 
constraints  with  planes;  contact  with  another  free  form  surface;  internal 
continuity  constraints  between  patches, 

• step  two  of  this  method  involves  mechanical  features.  These  features  are 
based  on  parameters  of  the  mechanical  model  (topology,  mobility  and 
force  density)  used  to  obtain  a deformed  geometry  and  on  the  choice  of 
a minimization  criterion.  Thus,  the  user  can  obtained  different  solutions 
with  a unique  set  of  geometric  features. 

These  first  two  steps  are  devoted  to  the  initialization  of  the  process.  The  user 
can  modify  one  or  all  these  features  if  he/she  does  not  accept  the  deformed 
geometry.  The  next  two  steps  are  transparent  for  the  user,  and  focus  on  the 
computation  of  the  deformed  geometry: 

• step  three:  a relationship  between  geometric  and  mechanical  features 
contributes  to  the  computation  of  new  external  forces  through  an  opti- 
mization process.  Different  minimization  functionals  can  be  incorporated 
into  this  process, 

• step  four:  these  new  forces  influence  the  static  equilibrium  positions  of 
bar  networks  (mechanical  models).  New  3D  positions  of  the  nodes  of  the 
bar  networks  are  computed,  i.e.,  new  3D  positions  of  the  vertices  of  the 
control  polyhedrons  of  the  trimmed  surfaces. 
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(a)  An  example  of  bar  network 


(b)  The  control  polyhedrons  of  Cal  and  Ca2  can  be  merged  to 
form  onebar  network  because  they  are  connected  along  their 


entire  boundary  whereas  Ca3  generates  another  bar  network 


since  Ca3  is  connected  to  Cal  and  Ca2  along  a trimming  line. 


Fig.  2.  Bar  networks  used  to  control  the  shape  of  surfaces. 


§3.  Mechanical  Deformation  Features 
Mechanical  model  of  deformation 

The  bar  network  (Figure  2a)  is  built  from  bars  B with  pin  joints  which  are 
assumed  to  rotate  without  friction  [4,11].  All  bars  are  under  tension.  The 
parameters  governing  the  static  equilibrium  state  of  such  a network  are 

• the  mobility  of  the  nodes  N,  i.e.,  fixed  (Nj)  or  free  (A))  to  move  in  3D 
space, 

• the  topology  of  the  bar  network,  i.e.,  the  way  the  bars  B are  connected  to 
the  nodes  N of  the  network, 

• the  force  density  qj  attached  to  each  bar  Bj  of  the  network  is  defined 
as  the  ratio  between  the  internal  force  fbj  into  the  bar  and  its  length  lj 

(qj  = qj  > 0).  The  positivity  constraint  ensures  the  tension  state 

in  every  bar  Bj , 

• the  external  force  f,  which  may  be  applied  to  the  ith  mobile  node  of  the 
bar  network. 


Linear  static  equilibrium 

Static  equilibrium  of  a bar  network  is  achieved  when  the  sum  of  the  external 
force  f,  applied  at  the  ith  node  equilibrate  the  forces  applied  by  each  bar 
meeting  at  that  node.  This  statement  becomes 

nbi 

fj  T ^ ] qj • (x/;  xn)  0,  i € {1, . . . ,tij],  (1) 

3=0 

where  x;,  stands  for  the  coordinates  of  the  ith  free  node,  qjf'x.k  - xu)  denotes 
the  tension  force  into  the  jth  bar  meeting  at  node  i,  iifn  designates  the  number 
of  bars  meeting  at  node  i and  n;  the  number  of  free  nodes  of  the  bar  network. 
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' JDi  .xi  + E)f  .Xf  - fx, 

< lDl.yl+  E)f  .yf  =fy, 

K K>i  .zi  + E)f  .Zf  = fz. 

for  one  bar  network,  where  E>i  is  a positive  definite  matrix  which  ensures  the 
existence  and  unicity  of  either  equilibrium  positions  (x;,  y;,  z;)  or  external 
forces  (fa,,  fy,  fz). 

Generally,  the  objective  is  to  calculate  the  3D  position  of  free  nodes. 
These  new  positions  are  obtained  with  a linear  equation  system  through 
changes  of  variables  of  the  parametric  process,  i.e.,  the  external  forces  f,. 

Analogy  used  for  surface  deformation 

The  deformation  method  uses  an  analogy  between  the  control  polyhedron  of 
a surface  and  the  mechanical  equilibrium  position  of  a bar  network.  One  bar 
network  is  associated  with  one  or  several  surfaces  as  follows: 

• the  nodes  of  a bar  network  coincide  with  the  entire  set  of  vertices  of  the 
control  polyhedron  of  a free-form  surface  (either  trimmed  or  not), 

• the  C°  continuities  are  directly  incorporated  into  the  mechanical  model. 
In  this  case,  several  bar  networks  can  be  merged  together. 


§4.  Deformation  Criteria 

The  geometric  constraints  generated  by  the  designer  to  prescribe  dimensions 
combined  with  the  equilibrium  equations  of  the  bar  networks  form  the  global 
set  of  constraints 


G — G(/lx,  • • • , /ra,:c,  /ly,  • • • , fniyi  flzt  • • • i fntz) i 

= G(F)  = 0,  i € nc}, 


(2) 


where  each  constraint  is  expressed  in  terms  of  the  external  forces  applied  to 
the  bar  networks.  External  forces  have  been  chosen  as  unknowns  rather  than 
force  densities  qj  because  they  produced  intermediate  solutions  which  reflected 
a real  deformation  process  of  a surface  whereas  iterating  with  qj  produced 
oscillations  around  the  solution.  Then,  the  resolution  has  been  conducted 
using  an  augmented  Lagrangian  method  to  provide  robustness  to  the  solving 
process. 

Because  the  number  nc  of  constraints  is  usually  significantly  smaller  than 
the  number  of  unknowns,  and  assuming  that  there  is  no  local  configuration 
with  an  overconstrained  subset  of  equations,  a functional  $(F)  can  be  associ- 
ated with  (2)  to  obtain  a solution  which  matches  a specific  designer’s  interest. 
Overconstrained  subsets  of  equations  can  be  encountered  when  a subset  of  (2) 
is  such  that  there  exists  locally  for  a given  surface  more  constraints  then  the 
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Fig.  3.  Relationships  between  external  forces  at  nodes  and  the  shape  of  a curve. 


number  of  free  parameters  xi,  yi,  zi.  When  G contains  nonlinear  equations, 
their  derivatives  are  related  to  the  geometrical  and  mechanical  parameters  by 

dGj  dGj  ^ ^ ^ ^ dG^  dZp  dzi^ 

13  dFj  dXp  ' dxlq  ' dFj  + dYv  ' dylq  ' dFj  + dZp  ' dzlq  ' dFj  ’ *' 

where  Jjp-,  are  related  to  the  geometric  constraints  set  by  the  de- 

QX  QY  dZ 

signer,  ggf-,  g^f~,  gjf  come  from  the  relationship  between  the  surface  and  the 
bar  networks  and  gfj-,  gjf,  gjfr  are  coefficients  of  Etf1. 

External  forces  at  nodes  and  shape  relationships 

According  to  (1),  external  forces  at  nodes  are  governed  by  the  length  of  the 
bars  as  well  as  the  angle  between  the  bars  meeting  at  a node.  Figure  3a 
illustrates  such  a configuration  for  a bar  network  which  corresponds  to  the 
control  polygon  of  a Bezier  curve  with  uniform  force  densities  in  its  bars. 

Then,  it  can  be  stated  that  a regular  bar  network  has  smaller  external 
forces  than  an  irregular  one,  since  regular  control  polygons  resemble  the  shape 
of  the  curve,  and  therefore  have  smaller  length  sides  and  wider  angles  between 
bars  (Fig.  3b). 

With  uniform  force  densities,  the  direction  of  the  external  forces  is  close 
to  the  bisecting  line  of  two  adjacent  bars  or,  for  a bar  network  attached  to 
a surface  patch,  this  direction  is  close  to  the  average  normal  direction  at  the 
given  node.  However,  this  behaviour  does  not  necessarily  generate  acceptable 
shapes  (as  depicted  in  Figure  4)  and  needs  to  be  combined  with  the  regularity 
criterion  of  the  intensities. 

A change  in  the  intensity  of  an  external  force  at  a node  N k results  in  a 
change  of  position  of  the  free  nodes  whose  amplitude  decrease  from  Nk  ■ The 
direction  of  movement  of  the  free  nodes  is  similar  to  that  of  the  external  force 
which  has  been  modified  [7]  at  Nk-  In  turn,  the  displacement  of  the  points 
on  a curve  or  on  a patch  follows  the  geometric  property  of  Bezier  or  B-Spline 
models,  i.e. , points  move  in  the  direction  of  the  movement  of  a control  point. 
Their  displacement  amplitude  is  therefore  smaller  than  that  of  Nk . 
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Fig.  4.  Specific  configuration  where  external  forces  coinciding  with  bisecting 
lines  does  not  provide  a smooth  curve. 


Fig.  5.  Criteria  related  to  the  external  force  at  a free  node. 

Various  deformation  criteria 

Based  on  the  previous  relationships,  three  categories  of  criteria  have  been 
identified: 

• the  first  one  is  related  to  the  external  forces  at  the  free  nodes  of  the  bar 
networks, 

• the  second  one  acts  simultaneously  over  the  external  forces  and  average 
normal  direction  at  the  free  nodes, 

• the  third  one  acts  over  the  average  force  around  a free  node,  but  is  not 
described  here  because  of  lack  of  space). 

The  first  category  takes  as  input  either  the  external  force  fjfc^  or  the  varia- 
tion Aff^of  this  force  at  the  iteration  [ k ] when  the  constraints  expressed  are 
nonlinear  (see  Figure  5). 

When  considering  the  functional  = min  ^TAF^  . AF^^ , the  de- 
signer expresses  the  minimum  change  in  the  shape  of  the  object  at  each  iter- 
ation until  the  constraints  G are  satisfied. 

When  considering  = min  (T(FW  - Fl0]) . (F[fc]  - F[01)j  as  the  func- 
tional associated  with  G,  the  designer  expresses  the  minimum  change  of  the 
object  shape  between  the  input  geometry  and  the  output.  In  case  of  linear 
constraints,  $[*],  and  are  identical.  This  functional  tends  to  preserve  as 
much  as  possible  the  previous  work  of  the  designer,  and  therefore  is  of  specific 
interest  during  a modeling  process. 
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Fig.  6.  Criteria  acting  over  the  intensity  and  the  direction  of  the  external  force 
at  a node. 


When  considering  the  functional  = min  ^rF^  . F^j , the  designer 
expresses  the  fact  that  the  resulting  surface  is  not  based  on  the  input  one 
since  the  initial  external  forces  are  not  taken  into  account.  Furthermore,  the 
minimization  of  the  intensity  of  the  external  forces  at  the  free  nodes  expresses 
that  the  output  control  polyhedrons  form  an  approximation  of  a minimal 
surface.  Hence,  the  output  surface  attached  to  the  bar  networks  represents  an 
approximation  of  the  minimum  surface  area  satisfying  the  constraints.  This 
approximation  is  even  more  effective  when  the  control  polyhedrons  converge 
toward  the  surface  itself. 


The  second  category  of  criteria  acts  simultaneously  over  the  intensity 
and  direction  of  the  external  forces  at  free  nodes  so  that  the  direction  of  the 
forces  meet  a given  criterion.  To  this  end,  an  average  normal  direction  is 


built  according  to  the  position  of  the  nodes  around  the  target  node.  Fig- 
ure 6 illustrates  the  planes  surrounding  the  ith  free  node  at  iteration  [A:  — 1] 
which  participate  to  the  definition  of  the  average  normal  direction 
When  considering  the  functional  i>pj0  = min  ^7  AF^0 . AF{fJ„  j , the  forces 
minimized  correspond  to  the  difference  between  the  projection  f^r~^  of  the 
external  force  f|fc  ^ onto  and  the  force  fj^  at  the  kth  iteration,  i.e. 


Ap[fcl  = ftfcI  _ 

pro— t i 


Using  this  criterion,  the  designer  expresses  that  the  output  surface  tends 


to  minimize  the  area  while  being  smooth  since  the  intensity  of  the  forces 
tend  to  decrease  like  as  well  as  the  direction  of  the  forces  tend  to  be  more 
regular  using  This  criterion  takes  into  account  the  shape  of  the  input 

geometry,  but  generates  a surface  which  is  ‘smoother’  than  with 

When  considering  the  functional  = min  ^TAF^t . AFJ.^(  j , the  forces 

minimized  correspond  to  the  difference  between  the  external  force  ro- 
tated onto  the  direction  of  and  the  force  ff'  at  the  kth  iteration, 

i.e.  AFj.^^  = fl^  — Ijfj*’  Such  a functional  tends  to  preserve 
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Initial  geometry 
Two  distinct  boundary 
conditions 


Initial  geometry 


Design  constraints 


Fig.  7.  Illustration  of  the  influence  of  the  deformation  criteria. 

the  intensity  of  the  external  forces  while  modifying  their  direction  in  order  to 
smooth  the  surface.  Combined  with  anisotropic  force  densities  in  the  bars  of 
the  networks,  this  criterion  expresses  the  deformation  behaviour  of  pipe-type 
objects  when  their  are  subjected  to  bending  deformations.  Though  the  previ- 
ous criteria  are  nonlinear  even  if  G is  linear,  their  efficiency  is  strong  enough 
to  justify  their  use  during  a design  process. 

§5.  Results  and  Examples 

The  above  deformation  criteria  have  been  applied  to  different  categories  of 
surfaces  to  illustrate  their  typical  behaviour  according  to  the  configurations 
described  in  the  previous  section.  Figure  7a  shows  the  effect  of  4>ioc,  7b  illus- 
trates 7c  $tot,  7d  $pro  and  7e  4>rot  under  various  designer  constraints. 

Two  distinct  input  geometries  are  used.  The  upper  one  is  a one  patch 
surface,  and  two  types  of  boundary  conditions  were  used,  i.e.  fixed  boundary 
lines  and  fixed  corner  points.  The  design  constraint  is  formed  by  a position 
and  tangency  constraint.  The  bottom  one  is  a multipatch  surface  with  G 1 
continuity  constraints,  where  the  designer  has  specified  position  and  tangency 
constraints  along  the  extreme  boundary  line  of  the  surface. 

§6.  Concluding  Remarks 

The  deformation  criteria  presented  provide  a diversity  of  control  of  shape. 
They  form  an  efficient  complement  to  the  geometric  constraints  set  by  the  de- 
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signer  to  let  him/her  adapt  the  result  to  his/her  needs.  Such  a diversity  cannot 
be  achieved  using  a mechanical  approach  solely  based  on  a minimization  of  the 
strain  energy  of  membrane  type  structures.  Furthermore,  the  criteria  set  up 
are  not  bound  by  a small  displacements  hypothesis  and  can  handle  geometric 
constraints  involving  significant  shape  changes. 

Future  work  will  focus  on  the  cross  influence  between  the  deformation 
criteria  and  the  boundary  conditions  applied  to  the  bar  networks  to  provide 
more  intuitive  user  interactions. 
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